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One can always decompose Dirichlet-Voronoi polytopes of lat- 
tices non-trivially into a Minkowski sum of Dirichlet-Voronoi 
polytopes of rigid lattices. In this report we show how one 
can enumerate all rigid positive semidefinite quadratic forms 
(and thereby rigid lattices) of a given dimension d. By this 
method we found all rigid positive semidefinite quadratic forms 
for d = 5 confirming the list of 7 rigid lattices by Baranovskii 
and Grishukhin. Furthermore, we found out that for d < 5 the 
adjacency graph of primitive L-type domains is an infinite tree 
on which GL(i(Z) acts. On the other hand, we demonstrate 
that in d = 6 we face a combinatorial explosion. 

1. Introduction 

Let L be a lattice in Euclidean space (M'^, (•,•))• With L one associates the 
Dirichlet- Voronoi polytope 

DV(L) = {x£R'^ : (x, x) < (v - x, v - x) for all v G L}. 

One can always decompose Dirichlet-Voronoi polytopes of lattices non-trivially 
into a Minkowski sum of Dirichlet-Voronoi polytopes of so-called rigid lattices. 
In this paper we show how one can find all rigid lattices in a given dimension. 

It will be more convenient to use the language of quadratic forms. With 
L we associate a positive definite quadratic form: By choosing a lattice basis 
bi, . . . , we get a positive definite matrix Q = {{hi, hj))ij giving the positive 
definite quadratic form x i— > x^Qx. 

By applying basic facts of Voronoi's theory of L-type domains we get an 
algorithm for finding all rigid forms of a given dimension. All quadratic forms 
whose Dirichlet-Voronoi polytopes have the same combinatorial-metric struc- 
ture belong to a pointed polyhedral cone, a so-called L-type domain. L-types 
domains of maximal dimension are called primitive, L-type domains of min- 
imal dimension 1 are called rigid, and elements of rigid L-type domains are 
called rigid forms. From Voronoi's algorithm for finding all primitive L-type 
domains we get the facets of every primitive L-type domain. By converting 
the facet description we find all extreme rays. If an L-type domain is invariant 
under a non-trivial symmetry group we can speed up the task of converting 
considerably using the adjacency decomposition method. 
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We applied this algorithm for dimensions < 5. Thereby, we confirmed the 
results by Baranovskii and Grishukhin (2): There is exactly one 1-dimensional 
rigid form (the Dirichlet-Voronoi polytope is a line segment), there are no 
2- and 3-dimensional rigid forms, there is exactly one 4-dimensional rigid 
form (the Dirichlet-Voronoi polytope is the 24-cell), and there are seven 5- 
dimensional rigid forms. Using this classification we verified that for d < 5 
the graph of primitive L-type domains is an infinite tree on which the group 
GLrf(Z) acts. It is computationally quite simple to perform these classifications 
because face lattices of these primitive L-type domains are very tame. 

This is no longer the case in dimension 6. We explored two notable primi- 
tive L-type domains of 6-dimensional quadratic forms. The first cone has 130 
facets and we do not know a primitive L-type domain in dimension 6 having 
more facets. The second cone has 100 facets and it contains a positive definite 
quadratic form associated to the lattice whose covering density is conjectured 
to be optimal in dimension 6. The automorphism groups of both cones are 
fairly big. By using the adjacency decomposition method we succeeded to 
compute the extreme rays of both cones: The fist one has 7, 145, 429 extreme 
rays, and the second one has 2, 257, 616 extreme rays. But many of these 
extreme rays correspond to equivalent (under the group GL^(Z)) rigid forms. 
In total we found 25, 263 non-equivalent rigid 6-dimensional positive definite 
quadratic forms. 

This paper is organized as follows: In Section [21 we recall some definitions 
and facts of Voronoi's theory of L-type domains. In Section |31 we show in 
which sense rigid forms are building blocks of Dirichlet-Voronoi polytopes. In 
Section we explain the adjacency decomposition method. In Section [S] and 
Section|21we report on computational results for rigid forms up to dimension 6. 

2. Notation: L-Type Domains 

In this section we define L-type domains for positive semidefinite quadratic 
forms. This enables us to define rigid forms. 

By S"^ we denote the space of all quadratic forms in d variables, and by 
<S>Q we denote the set of all positive semidefinite quadratic forms which is a 
closed pointed cone. In the following we will identify S'^ with the space of 
all symmetric (d x (i)-matrices. We say that two quadratic forms Q,Q' are 
arithmetically equivalent if there exists an integral unimodular matrix A € 
GLrf(Z) so that Q' = A*QA. 

Let Q G 5>Q be a positive semidefinite quadratic form arithmetically equiv- 
alent to ( Q ) where Q' is positive definite. We define the Dirichlet- Voronoi 
polytope of Q by 

DV(Q) = {x*Q G {W^y : x*Qx < (x - v)*Q(x - v) for ah v G Z'^]. 

This way of defining Dirichlet-Voronoi polytopes of positive semidefinite 
quadratic forms is due to Namikawa (7). The definition has the important 
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feature that we can define L-type domains of positive semidefinite quadratic 
forms by using the concept of strongly isomorphic polytopes which we recall 
now. 

Let y be a finite dimensional real vector space. Let P C 1/ be a convex poly- 
tope, and let / G be a linear functional. We define the support functional 
of P by r?(P, /) = max{/(a;) : x G P} and by P/ = {x G P : f{x) = r]{P, /)} 
we denote the face of P in direction /. Wc say that two convex polytopes 
P,Q & V are strongly isomorphic if for every f,g € V* with Pf C Pg we have 

QfQQg. 

We say that two positive semidefinite quadratic forms Q, Q' belong to the 
same L-type if their Dirichlet-Voronoi polytopes are strongly isomorphic. The 
set of all positive semidefinite quadratic forms belonging to the same L-type 
is called an L-type domain. Two L-type domains A, A' are arithmetically 
equivalent if there exists A G GLrf(Z) so that A' = A^AA. 

In (11) Voronoi showed that L-type domains are open pointed polyhedral 
cones, that the L-type domains give a face-to-face partition of 5>q, and that 
there are only finitely many non-equivalent L-type domains. L-type domains 
which arc of maximal dimension ^^'^^^^ arc called primitive. L-typc domains 
which are of minimal dimension 1 are called rigid. Positive semidefinite qua- 
dratic forms lying in a rigid L-type domain are called rigid, too. 

3. Rigid Forms and Dirichlet-Voronoi Polytopes 

In this section we show in what sense Dirichlet-Voronoi polytopes of rigid 
positive semidefinite quadratic forms are building blocks of Dirichlet-Voronoi 
polytopes of general positive semidefinite quadratic forms. 

Lemma 1. Every Dirichlet-Voronoi polytope of a positive semidefinite qua- 
dratic form is Minkowski sum of Dirichlet- Voronoi polytopes of rigid forms. 
More precisely: Let A be the topological closure of an L-type domain. For pos- 
itive semidefinite quadratic forms Qi,. ■ ■ ,Qn £ A and non-negative numbers 
ai,. . . ,an we have 

n n 

DV(^aigi) = I]«iDV(Qi). 

i=l i=l 

The authors do not know exactly the origin of this lemma. Loesch gave it 
in a dual formulation in (6). Later, Ryshkov gave in (9) a similar but less 
precise statement. 

4. Computational techniques 

Fukuda's program odd (5) computes the list of extreme rays of a polyhedral 
cone given its list of facets. In our case the number of extreme rays can be very 
large so that we cannot use odd naively. We apply another technique called 
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adjacency decomposition method to use the symmetry of the polyhedral cones 
we are considering. 

Let C C R*^ be d-dimcnsional polyhedral cone determined by a set of facets 
We assume that C is pointed at the origin. By fi G {W^)* 
we denote a linear functional defining Fj, i.e. Fj C {x G M'' : /i(x) = 0}. 
Let E be an initial extreme ray of C which we find e.g. by solving a generic 
linear program on C. Wc compute the extreme rays adjacent to E: first we 
project C along F by a linear map tt. Then, we find the extreme rays of 
this projected cone 7r(C) (using cdd, or applying this procedure recursively). 
Every extreme ray F^ of 7r(C) corresponds to a two-dimensional face F of C 
in which F lies. Therefore, there is exactly one more extreme ray E' of C in 
F. Every e' € F' can we written as e' = ae + /Jbtt, with e G F, G F^ and 
some a,/3, which we have to compute. This can be done by solving a "two- 
dimensional linear program": fi{ae + /Se.,^) = for all facets of C incident to 
F, fj{ae + /^Btt) > for all other facets of C. The key computational step 
in the procedure above is the computation of the extreme rays of 7r(C). The 
complexity of this computation is related to the incidence number of F, i.e. 
the number of facets containing F. 

The adjacency decomposition method applies to polyhedral cones having a 
non-trivial symmetry group: 

(1) Take an initial list of orbits of extreme ray of C. 

(2) Take a representative F of an orbit and finds the extreme rays (Fj)i<j<TO 
adjacent to it. 

(3) If some Fj represents a new orbit, then we add it to the list of orbits. 

(4) Finish when all orbits have been treated. 

This procedure has two main computational bottlenecks: it can be difficult to 
identify new orbits, and the incidence of extreme rays can be too high. 

Since the symmetry groups of the cones we considered was not too big, the 
first bottleneck was not a problem: for every new extreme ray we generated 
the whole orbit. 

For dealing with the second bottleneck we used Balinski's theorem: 
Theorem 1. ((1), see e.g. (12)) 

Let C he a d- dimensional pointed polyhedral cone. Let G he the undirected 
graph whose vertices are the extreme rays of C and whose edges are the 2- 
dimensional faces of C . Two vertices Fi,F2 are connected by an edge F if 
Fi, F2 G F. Then, the graph G is {d— l)-connected, i.e. removal of any d — 2 
vertices leaves it connected. 

Due to Balinski's theorem, we can replace the criterion "Finish when all 
orbits have been treated" by "Finish when the number of extreme rays in 
untreated orbits is lower than d — 2" . 

Wc applied the adjacency decomposition technique, starting with orbits of 
lowest incidence. After some time, we found the complete list of orbits. But 
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still we had to treat the orbits with highest incidence. These very degenerate 
orbits have usually a particular significance, in our case they correspond to 
quadratic forms lying in the boundary of 5>o- In the L-type domain con- 
sidered, we find out a peculiarity: the number of elements of the orbits with 
highest incidence is very low, actually lower than d—2. Therefore, those orbits 
cannot disconnect the skeleton graph and so we can stop earlier avoiding the 
computation of adjacencies of those orbits. 

5. Dimensions 1, . . . , 5. 

In dimension 1 all positive definite quadratic forms are rigid; their Dirichlet- 
Voronoi polytopes are line segments. In dimension 2 and 3 there are no rigid 
positive definite quadratic forms because in these dimensions every Dirichlet- 
Voronoi polytope is a zonotope whence it is a Minkowski sum of line segments. 
In all these dimensions there is only one non-equivalent primitive L-type do- 
main. In dimension 4 there are three non-equivalent primitive L-type domains, 
and there is exactly one rigid positive definite quadratic form which is asso- 
ciated to the root lattice D4. Its Dirichlet-Voronoi polytope is the 24-cell. In 
dimension 5 there are 7 rigid positive definite quadratic forms. First, they 
were enumerated by Baranovskii and Grishukhin (2) by using Engel's list of 
zone-contracted lattices (4). Our computations confirmed their result. 

Now we will argue that these computations do not require much computa- 
tional effort. All L-type domains up to dimension 4 are simplicial polyhedral 
cones. As noticed first by Barnes and Trennery (3) (see also the discussion 
in (8) §13) this docs no longer hold in dimension 5 and above. The follow- 
ing table shows how the numbers of facets of primitive L-type domains are 
distributed among the 222 non-equivalent primitive L-types domains in di- 
mension 5. With n we denote the number of facets and with Li(n) we denote 
the number of non-equivalent primitive L-type domains in dimension 5 having 
exactly n facets. 
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The next table shows the distribution of the numbers of extreme rays among 
the 222 non-equivalent primitive L-type domains in dimension 5. With n we 
denote the number of extreme rays and with L2(n) we denote the number of 
non-equivalent primitive L-type domains having exactly n extreme rays. 
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The following table shows how the ranks of the extreme ray are distributed 
among the 222 non-equivalent L-type domains in dimension 5. By n we denote 
the number of extreme rays of and by Rk{n) we denote the number of non- 
equivalent L-type domains having exactly n extreme rays containing positive 
semidefinite quadratic forms of rank k, k € {1,4,5}. 
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By our computations we found out that for d < 5 the graph of primitive 

L-typc domains (vertices = primitive L-type domains, edges = facets between 
primitive L-type domains) is an infinite tree on which the group GL(i(Z) acts. 
The grapli of primitive L-type domains is a tree if and only if it has no cycle 
and therefore if and only if every ridge contains at least one degenerate form. 
Is this always the case? We think that the answer is "No" even for d = 6. We 
also think that there is a primitive L-type domain whose extreme rays are all 
non-degenerate. 

6. Dimension 6 

We consider here two different primitive L-type of 6-dimensional positive 
semidefinite quadratic forms, which were considered by the second author in 
(10). 

6.L The Cone Ci. The cone Ci is a primitive L-type domain of 6-dimensional 

positive semidefinite quadratic forms. It has 130 facets and we conjecture 
that there is no L-type domain of 6-dimensional positive semidefinite qua- 
dratic forms having more facets. The automorphism group of Ci has order 
1920. Using the adjacency decomposition method we computed that Ci has 
7, 145, 429 extreme rays in 4, 440 orbits. It is worthwhile to note that among 
the 4, 440 orbits there are two orbits which are equivalent under the group 
GLq{Z), and that there is one orbit which does contain forms of rank 5, so 
that we found altogether 4, 438 non-equivalent (under the group GL6(Z)) rigid 
positive definite quadratic forms. 

6.2. The Cone C2. The cone C2 contains the 6-dimensional positive definite 
quadratic form associated to the best known 6-dimensional lattice covering 
(10). It has 100 facets and its automorphism group has order 120. Using the 
adjacency decomposition method we computed that C2 has 2, 257, 616 extreme 
rays in 20, 871 orbits which correspond to 20, 861 non-equivalent (under the 
group GLq{Z)) rigid positive definite quadratic forms. Three orbits correspond 
to positive semidefinite quadratic forms being not positive definite. 

6.3. Connection between Ci and C2. In this section we show how Ci and 
C2 are related. Both cones Ci and C2 contain the rigid form 
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associated to the lattice Eg. The automorphism groups of Gj of Cj, z = 1, 2 are 
subgroups of the automorphism group G = {T £ GL6(Z) : T^Q^*T = Qe*} of 

The subspace Ii of all quadratic forms invariant under the group Gi is 
spanned by Qe* and Ri (see below) which is an extreme ray of Ci. If we 
intersect the cone G2 with Ii we get a two-dimensional cone with extreme 
rays Qe;; and i?2- The rigid forms Ri and R2 are 



Ri = 



(12 


3 


6 


6 


-3 


3\ 




/O 














0\ 


3 


7 


4 


4 


3 


2 







5 


2 


2 


3 


1 


6 


4 


8 


3 


1 


4 







2 


4 





2 


2 


6 


4 


3 


8 


1 


4 





2 





4 


2 


2 


-3 


3 


1 


1 


7 


3 







3 


2 


2 


5 


3 


V3 


2 


4 


4 


3 






^0 


1 


2 


2 


3 


5/ 



6.4. Further Remarks. It is remarkable that there are only very few in- 
stances of arithmetically equivalent extreme rays, which are not equivalent 
under the symmetry group of the cone. The number of extreme rays of both 
cones is extremely large. Nevertheless it is interesting to note that the number 
of non-equivalent extreme rays corresponding to forms which are not positive 
definite is low. In total, we obtained 25, 263 new rigid positive definite qua- 
dratic forms in dimension 6. 
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